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We analyze the magnetic excitations of a spin-1/2 antiferromagnetic Heisenberg 
model with alternating nearest neighbor interactions and uniform second neighbor 
interactions recently proposed to describe the spin-Peierls transition in CuGe03. 
We show that there is good agreement between the calculated excitation dispersion 
relation and the experimental one. We have also shown that this model reproduces 
satisfactorily the experimental results for the magnetization vs. magnetic field curve 
and its saturation value. The model proposed also reproduces qualitatively some 
features of the magnetic phase diagram of this compound and the overall behavior 
of the magnetic specific heat in the presence of applied magnetic fields. 
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Interest in low-dimensional quantum spin systems has been intensified in recent years by 
the observation of a rich variety of behaviors in new materials and by important theoretical 
results, In particular, the spin-Peierls (SP) transition is one of the most attractive 
phenomena and it was first discovered and studied in organic compounds. |3[] At the SP 
transition, the high-temperature phase with uniform spin-1/2 Heisenberg antiferromagnetic 
chains ("U phase") changes into the low-temperature phase with dimerized or alternating 
chains ("D phase"). In the dimerized phase the system is characterized by the presence of 
an energy gap (the "SP gap" ) between the singlet ground state and the first triplet state. 

Recently the SP transition was observed in the inorganic material CuGeO"3 by Hase, 
Terasaki and Ichinokura. [|J They showed that the magnetic susceptibility in CuGe03 along 
all directions drop drastically at the transition temperature T§p « 14 K. The SP nature of 
this transition has been confirmed by several experimental results. [^||6|,|7|,|8|,|9|,p!0[| The quasi 
one-dimensional (ID) nature of this compound is given by the presence of weakly coupled 
linear Cu 2+ spin-1/2 chains along the c direction. The transition is the result of the coupling 
between these chains and the three-dimensional phonons. ||[7| 

In order to characterize the SP phase in this compound, it is important to investigate 
its thermodynamic properties in the presence of magnetic fields because as it is well-known 
in organic SP materials they have characteristic behaviors below Tgp and large magnetic 



fields. [|TTJ The phase diagram of SP systems in the magnetic field and temperature plane, 
as determined experimentally from specific heat and ac susceptibility measurements, shows 
the uniform and dimerized phases at low magnetic fields, and another magnetic phase ( "M 
phase") at high magnetic fields. In the absence of magnetic fields, the dimerization of 
the chains corresponds to a lattice distortion with wave vector it. For any T < Tsp, the 
magnetic field above a critical value H C (T) favors a distortion at some other wave vector or 
another type of distortion. The M phase could then correspond to another conmensurate, 
discommensurate (or magnetic soliton) or inconmensurate phase. In the case of CuGe03, 
this M phase has not yet been fully characterized. 



Recently, in Ref. |12| ("RD" in the following) we attempted a description of the SP 



transition in CuGeO"3 with a spin-1/2 ID Heisenberg antiferromagnetic model coupled to 
the three-dimensional phonons in the adiabatic approximation. The spin part of this model 
is defined by the Hamiltonian: 

H B = X)J i S i .S i+1 + aJX;Si-S i+2 (1) 

i i 

where the index i runs over the lattice sites (i=l,...,N, N: number of sites) with periodic 
boundary conditions. The exchange coupling constant J; = J(l + S(— l) 1 ) where 5 = S(T) 
is a temperature dependent dimensionless quantity which gives the degree of dimerization. 
The model Eq. ([!]) has been studied so far in the limits a = fli~3 |, and 5 = 0. |L 



The parameters of the model were chosen so as to give a good fit of the experimental 
magnetic susceptibility and to reproduce the measured excitation gap at T = 0. 0,§] The 
resulting parameters are: J = 160 K, a = 0.36 and 8(T = 0) = 0.014. As it is well-known, 
the ID Heisenberg antiferromagnetic model with first and second neighbor interactions has 



a finite singlet-triplet gap for a > a c ~ 0.2412. fnjTS] Thus, the value of a = 0.36 implies 
that there is a finite gap even in the absence of dimerization. This spin gap was estimated 
in RD to be ~ 0.015 in units of J, or 2.4 K, much smaller than the smallest measured value 
for CuGe03. 0] In a more recent study Castilla, Chakravarty and Emery IfTSJ argued, from 
an analysis of the neutron scattering measurements of the superlattice peaks @, that the 
value of a should be below the critical one. For this reason, they proposed the following set 
of parameters: J = 150 K, a = 0.24, and 5(T = 0) = 0.030. However, we think that this 



argument, based on the theory of Cross and Fisher |16j which uses a Luther- Peschel-type 



treatment of the spin Hamiltonian, is not conclusive. Moreover, we should notice that there 
is a rather large interchain coupling along the b direction, Jb/J ~ 0.1 that could lead to 
a different effective a c . 

The purpose of this work is to describe using the model proposed in RD other features of 
the SP phase in CuGeOs, in particular the nature of magnetic excitations and the behavior 
of various thermodynamic properties in the presence of an external magnetic field as have 
been revealed by an increasing amount of experimental results. We want to show also that 



the model and numerical procedure developed in RD could provide an alternative complete 
theory of the spin-Peierls transition and that it could be considered as an useful tool to 
analyze experimental data. These studies will be performed by using exact diagonalization 
techniques on finite chains. In these calculations we adopt the parameters given in RD which 
lead to a somewhat better fit to the experimental magnetic susceptibility as well as for the 



dispersion relation as we shall see below than the parameter set given in Ref. [fLql . Most of 
the calculations were performed also for this latter parameter set and we shall discuss the 
main qualitative differences found. 

We start our study with an analysis of magnetic excitations at zero temperature and in 
the absence of a magnetic field. For this study we have employed the Lanczos algorithm for 
lattices of up to 24 sites. We have also calculated the dynamical magnetic structure factor 



S(k, u) for all the wave vectors k. [17 This dynamical response has the largest intensity at 
k = (corresponding to ir in the uniform lattice) and hence we shall concentrate our study 
to this momentum. The results for S(k = 0,uj) are depicted in Fig. la for a = 0.36 
and S = 0.014. We also show the results for the nondimerized model for comparison. The 
first peak corresponds to the singlet-triplet excitation. This peak has the largest weight, in 
agreement with experimental results. || We also note other excitations at higher frequencies 
and with lower intensities. It is interesting to notice that the second peak in the nondimerized 
case is split by a finite 5. The relative intensity of this second peak with respect to the first 
one remains approximately constant as N is increased. Similar results are obtained with 



the set of parameters of Ref. |[15|| , although the second peak is located at slightly higher 
frequencies than for our set of parameters. In order to understand the origin of these 
excitations, let us examine them as a function of k. The results for the N = 24 chain are 
shown in Fig. lb together with recent neutron scattering results. || It can be seen that the 
parameter set of RD fits slightly better the experimental data than the parameter set of Ref. 
fT5|| . In Fig. lc we show the excitation energy of various states as a function of k for the 
same chain. At k = it is convenient to label each energy level by E ; (S, P) where S is the 
total spin, P (= +1, —1) is the quantum number associated to the reflection symmetry, and 



i = 0(1) for the ground state (first excited state). For N = 4n, n integer, the singlet-triplet 
spin gap is A st = E (l, — ) — E (0, +). For N = 4n + 2, the values of P are reversed. The 
values of P quoted below correspond to the case N = 4n. 

By studying the excitation spectra for several lattice sizes, we obtained that the position 
of the second peak in S(k = 0,uj) is uj\ = Ei(l, — ) — E (0, +). The excited singlet state 
Ei(0, +) appears inside the singlet-triplet gap but this is just a finite size effect as can be 
seen in Fig.ld. This feature is absent for a = 0.24, i.e. Ei(0, +) > E (l, — ) for any lattice 
size. In Fig. Id, we plot A st , uj\ and u 2 = Ei(0, +) — E (0, +) as a function of N _1 . The 
data do not seem to vary linearly with N -2 and there is some spread in the extrapolated 
values depending on the fitting function chosen. In Fig. Id tentative extrapolations using 
the function u = + aN -1 + bN -2 are shown. It is apparent that in the bulk limit the 
singlet-triplet excitation becomes the lowest excitation in agreement with the case a = 0.24 
and with experimental results. Besides, several trial extrapolations indicate that there would 
remain a finite gap between the first and second peaks of S(k = 0, to) in the thermodynamic 
limit for a = 0.36. The results are less conclusive for a = 0.24. Finally, we want to notice 
that the first excited triplet state has P = — 1 for a = 0.36 but P = +1 for a = 0.24. 

We turn now to the study of thermodynamic properties in the presence of an applied 
magnetic field H along the c direction. In order to study the SP phase at finite temperature 
and in the presence of a magnetic field we minimize the free energy for the total Hamiltonian 
with respect to 5 to obtain 5 eq (H, T). The total Hamiltonian of the system consists of the 
spin part given by Eq. (|l|) with the addition of an elastic term ~ 5 2 and the Zeeman term 
g/! S z H, where g is the Lande factor (we take g = 2), /x is the Bohr magneton and S z is 
the z-component of the total spin. The dimerization parameter 8 eq (H, T) is, for a given T, a 
decreasing function of H (Fig. 2a). Consequently, the spin gap and the critical temperature 
T SP (H) are also suppressed by the application of an external magnetic field, as discussed 
in several theoretical studies. [|11].[1(^.|T^] Once <5 e(? (H, T) is known, the computation of the 
thermodynamic properties we are interested in is straightforward. This procedure only 
requires the computation of the eigenvalues of the Hamiltonian ([!]) which is readily done by 
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using the Householder algorithm for chains of up to 14 sites. 

Recently the magnetization curve was measured in the presence of ultra-high magnetic 



fields up to ~ 500 T. ||T9[ The measurement of the magnetization and its saturation can be 
used to determine the values of the antiferromagnetic exchange interactions for an adopted 
model Hamiltonian. Nojiri et al Jill noticed that the experimental magnetization curve 



shows good agreement with the theoretical curve obtained for the spin-1/2 uniform Heisen- 
berg antiferromagnetic chain with nearest-neighbor coupling J = 183 K. However, as dis- 
cussed in RD, this model does not satisfactorily reproduce the magnetic susceptibility data. 
We plot in Fig. 2b the magnetization M versus applied magnetic field H obtained with 
the model proposed in RD for N = 14 site chain at T = 5 and 8 K. We also include the 



experimental results corresponding to T = 6 and 10 K. |19] Our numerical curves present 
the typical steplike structure due to the finite lattices involved. We observe an overall 
good agreement of the theoretical curve compared to the experiment results. The temper- 
ature dependence of our results is somewhat weaker than the experimental one, specially 
for high magnetic fields where presumably the three-dimensional effects are more impor- 
tant. We point out here that without further parameters to adjust in our model, we have 
also obtained the value of the saturation field (~ 253 T) in reasonable agrement with the 
experimental result. We obtain very similar results with the parameters of Ref. fl5| . 



The knowledge of 5 eq (H, T) enables us to discuss some features of the magnetic phase 
diagram of CuGeC>3. This magnetic phase diagram has been partially determined experi- 
mentally by measuring the magnetization in magnetic fields up to 25 T, || and it presents 
the three phases characteristic of a SP system. Hase et al. || showed that as the applied 
magnetic field is increased, at temperatures below Tsp(0), there is a nonlinear increase of the 
magnetization above a certain H C (T). |22| This nonlinear increase, together with the pres- 
ence of hysteresis, indicates that there is a first order transition between the dimerized phase 
and the high-field M phase for T < 10 K. The hysteresis disappears above this temperature 
and the D-U phase boundary becomes a second order one. This boundary is determined in 
our microscopic theory by the condition 5 eq (H, T) = 0. The results for N = 12, 14 and 16 



site chains [^T] using our parameters set are plotted in Fig. 3. We also show for comparison 



the analytical results obtained by Cross |23[ using the theory developed in Ref. |lb| . The 
most important assumption of this theory is that Tgp(O) <c J, moderately satisfied in this 
compound. For low magnetic fields, our results closely agree with the curve calculated by 
Cross corresponding to the D-U transition (dashed line) and with experimental results. PJ2U|] 

For higher magnetic fields, the theory of Ref. predicts a transition between the high- 
field magnetic (assumed incommensurate) and the uniform phases, which is also shown in 
Fig. 3 (full line). There are still no experimental results for CuGe03 corresponding to this 
M-U transition. In order to study the incommensurability of the M phase and the transitions 
D-M and M-U, the nearest neighbor coupling constants in our model should be generalized 
as Ji = J(l + <5cos(qi)) where q = q(H) is the component of the phonon wave vector parallel 
to the chain direction. The dimerization studied so far corresponds to the commensurate 
value q = 7r. In the presence of a magnetic field the value of q is the one that minimizes 
the free energy of the system. At present there are important numerical limitations to 
accomplish such a program. Moreover, it should be noticed that the interchain exchange 
interactions have stronger effects in the M phase. [[24]] Although we have not attempted an 
extrapolation to the bulk limit, the curves obtained with our model in the high-field region 
should converge as the chain size is increased to an upper bound to the D-M boundary. The 
determination of the first-order nature of the D-M phase transition is also out of the scope 
of this work because of the limitation in the numerically accessible chain sizes. 

Finally, we study the magnetic specific heat C m in the presence of a magnetic field. The 
results for C m /T obtained for the N = 14 site chain as a function of T/Tgp for H = 0, 5 
and 10 T are shown in Fig. 4. Experimentally, it has been found that the specific heat 



shows a sharp anomaly at Tsp due to the SP transition p5|j2q] . This peak is clearly seen 
in our results. Consistently with the results shown in Fig. 2b, the position of the peak 
corresponding to the SP transition decreases as the magnetic field is increased and the low 
temperature part of the curve moves upward. These changes of C m /T with magnetic field 



are also in qualitative agreement with the experimental data from Ref. |26| . The differences 



in shape between the theoretical and experimental curves are presumably also a finite size 
effect. 

In summary, we have analyzed magnetic excitations of the spin-1/2 antiferromagnetic 
Heisenberg model with alternating nearest neighbor interactions and uniform second neigh- 
bor interactions proposed in RD to describe the spin-Peierls transition in CuGe0 3 . These 
results are among the first reported for this very complicated model and they could be 
relevant for future neutron scattering experiments. There is good agreement between the 
calculated excitation dispersion relation and the experimental one. We have also shown that 
this model reproduces satisfactorily the experimental results for the magnetization vs. mag- 
netic field curve and its saturation value. The model proposed also reproduces qualitatively 
some features of the magnetic phase diagram of this compound and the overall behavior of 
the magnetic specific heat in the presence of applied magnetic fields. 

We acknowledge H. Nojiri for sending us their results for M vs. H prior publication. S. 
K. acknowledge helpful discussions with A. Dobry and A. Greco during the early stages of 
this work. 
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FIGURES 

FIG. 1. a) Spin dynamical structure factor S(k = 0,to) for N = 20 with a = 0.36. The solid 
(dashed) line corresponds to 5 = 0.014 (5 = 0). b) Excitation dispersion for the N = 24 chain. The 
solid line corresponds to J = 160 K, a = 0.36 and 5 = 0.014 and the dashed line to J = 150 K, 
a = 0.24 and 5 = 0.030. The experimental data (open circles) is taken from Ref. [g]. c) Excitation 
dispersion of various states for the N = 24 chain with J = 160 K, a = 0.36 and 5 = 0.014. The 
two lowest singlet (triplet) states are indicated with squares (circles), d) A s t (solid circles), uj\ 
(squares) and UJ2 (open circles), defined in the text, versus N _1 . The dotted lines are quadratic fits 
in N- 1 . 

FIG. 2. a) The reduced dimerization constant 5(H) versus the applied magnetic field H at 
T = 3 K (circles), 5 K (squares), 7 K (diamonds) and 9 K (triangles), b) The magnetization M 
as a function of magnetic field at different temperatures. The theoretical results for N = 14 are 



plotted with lines and the experimental data from Ref. [19| with symbols. 



FIG. 3. Magnetic phase diagram of the SP system CuGeOs. Open circles, solid squares and 
open diamonds correspond to S eq (T, H) = for N=12, 14 and 16 sites respectively. We also show 
the experimental data from Ref. [|| (crosses). The results obtained by Cross [^| are plotted with 
solid and dashed lines. 

FIG. 4. Magnetic specific heat as a function of temperature, for different values of the applied 
magnetic field H. 
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